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Nonequilibrium steady states in an open system connecting two reservoirs of platelike colloidal 
particles are investigated by means of a recently proposed phenomenological dynamic density func- 
tional theory [M. Bier and R. van Roij, Phys. Rev. E 76, 021405 (2007)]. The plateUke colloidal 
particles are approximated within the Zwanzig model of restricted orientations, which exhibits an 
isotropic-nematic bulk phase transition. Inhomogeneities of the local chemical potential generate a 
diffusion current which relaxes to a nonvanishing value if the two reservoirs coupled to the system 
sustain different chemical potentials. The relaxation process of initial states towards the steady 
state turns out to comprise two regimes: a smoothening of initial steplike structures followed by an 
ultimate relaxation of the slowest diffusive mode. The position of a nonequilibrium interface and 
the particle current of steady states depend nontrivially on the structure of the reservoirs due to 
the coupling between translational and orientational degrees of freedom of the fluid. 
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; I. INTRODUCTION 

Complex fluids of platelike colloidal particles, e.g., 
clay suspensions, have been investigated to quite an 
extent in recent years. The scientific interest largely 
stems from the enormous range of phenomena found 
in these systems such as flocculation, glass transitions, 
gelationjaging, and even li quid crystal phase transitions 

■ Ens mi la S S 1 MM, M due to the inter- 
play of translational and orientational degrees of freedom 
of the constituting particles. Whereas the liquid crystal 
properties of fluids of platelike colloidal particles have 

' been addressed in some theoretical studies devoted to 
homogeneous as well as inhomogeneous equilibrium sys - 

tems 14 M M [13 [H, [H, H mi HI i2a, m, M m , 

not many investigations of the nonequilibrium behavior 
have been performed. 

In a recent publication the authors proposed a 
phenomenological dynamic density functional theory 
(DDFT) in order to describe the relaxation dynamics 
of fluids of platelike colloidal particles under the influ- 
ence of an external field [l^l- Relaxation into stable 
' or metastable states, both characterized by a homoge- 
1 neous local chemical potential, has been found depend- 
ing on the initial state and the external field. In the 
present contribution, applying the same DDFT formal- 
ism, nonequilibrium steady states are investigated which 
form within a channel connecting two particle reservoirs 
of different chemical potentials. As the local chemical 
potential in the channel is expected to be spatially inho- 
mogeneous DDFT can be considered as the natural for- 
malism to describe the fluid of platelike colloidal particles 
on a coarsegrained level. 
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DDFT is an extension of equilibrium density functional 
theory (DFT) [H, [2^, Is^l to nonequilibrium conditions 
by propo sing an equation of motion for the one-particle 
densities [31|. The one-particle densities are assumed to 
describe the (dynamical) state of the system completely. 
On the one hand, the DDFT equations are similar to the 
traditional time-dependent Landau- Ginzburg and Cahn- 
Hillard models of critical dynamics , sp inodal decompo- 
sition, and crystal growth [H, IHj IHj 113 ■ 
the other hand, the DDFT equations have been derived 
within the framework of (overdamped) Langevin dynam- 
ics, which is considered a reasonable description for di- 
lute colloidal dispersions [H, [s^ HO, Hlj . As the number 
densities of the present work are close to the isotropic- 
nematic two-phase region, which is located at small den- 
sities for highly anisotropic particles, the applied DDFT, 
which neglects hydrodynamic interactions, is expected to 
be valid. The DDFT proposed in Ref. [27] is briefly sum- 
marized in Sec. [ill 

The present work is restricted to fluids of platelike col- 
loidal particles, i.e., nonergodic states such as glasses or 
gels are beyond the scope of this contribution. It is known 
that an ergodic system described by a linear master equa- 
tion relaxes towards a unique steady state [i^l . Although 
the formalism applied here is not equivalent to a linear 
master equation, the relaxation towards a unique steady 
state is nonetheless expected. This relaxation process is 
described in Sec. Illll in terms of the number density, the 
orientational order, the local chemical potential and the 
particle current. 

The final steady state is necessarily a nonequilibrium 
steady state, because the particle reservoirs coupled to 
the system preclude equilibration of the system. This 
nonequilibrium steady state exhibits a nonvanishing par- 
ticle current which is sustained by the chemical potential 
difference of the particle reservoirs coupled to the system. 
Moreover, for suitably chosen reservoir chemical poten- 
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tials, the nonequilibrium steady state shows signs of bulk 
phase transitions. These issues are discussed in Sec. IIVI 
Section |V] discusses the results found in Sees. |TTT] and 
IIVI and closes with a short summary. 



II. FORMALISM 
A. Model fluid and system geometry 

Consider a dispersion of monodisperse, hard, infinitely 
thin, square colloidal particles within a three-dimensional 
continuous solvent. The side length of the square parti- 
cles is D. The orientations, described by the normal vec- 
tor of the square face, are restricted to directions parallel 
to the Cartesian axes (Zwanzig approximation |43|). A 
particle is called an i particle if its orientation is along 
the i axis, i £ {x, y, z}. 

The system under consideration is a channel of length 
H which connects two particle reservoirs. The channel is 
assumed to be much wider than the particle size D such 
that effects of the channel walls onto the colloidal fluid 
are negligible. Consequently the fluid structure in the 
channel is expected to vary only along the channel axis 
which is taken as the z axis. The channel is located at 
the z axis interval [0,-ff]. The local number density of i 
particles in the channel at position z e [0, i?] is denoted 
by gi{z). The abbreviation g :— {Qx, Qy, Qz) is used later. 

The structure of the model fluid is adequately de- 
scribed in terms of the total density g := gi and 
the order parameter tensor Q [4^ 1 which, within Zwanzig 
models, is given by 



Qii' — 2 ^ ] ' 



(1) 



where Sw is the Kronecker delta. 

The two particle reservoirs are connected to the chan- 
nel at positions z = and z — H, respectively. They 
are assumed to sustain equilibrium bulk structures g 
and g^ corresponding to the chemical potentials /iq a-nd 
fiH, respectively. The coupling of the two reservoirs to 
the system amounts to the Dirichlet-like boundary con- 
ditions g{z < 0,t) = g^ and g(z > H,t) — g^. 



B. Dynamic density functional theory 

Within this work the dynamic density functional the- 
ory (DDFT) proposed by the authors in Ref. [131 is ap- 
plied. It consists of the set of equations of motion for the 
one-particle densities gi given by 



dgi{z,t) 
dt 



/ dgi{z,t) 
\ dt 



/ dg.,{z,t) 
V dt 



(2) 



with the translational diffusion for fixed orientation de- 
scribed by 



/ dg^{z,t) 
\ dt 



dz 



where the particle currents along the z axis are 

d(3fii{z, [g]) 



[g]) -T,g,{z)- 



dz 



(3) 



(4) 



and with the rotational diffusion at fixed position z de- 
scribed by 



/ dgi{z,t) 
V dt 



(^(3n,{z,[g{t)]) - f3^i^,{z,[g{t)])] 



The driving force is due to inhomogcneitics of the local 
chemical potential 



l^t{z, [q\) 



6F 



Sgi{z) 



(6) 



which derives from the free energy functional [45l . lit 
(3F[£\ = /dz (^^ft(z) (ln(ft(z)A3) - l) + 

Hn{z))' 



(7) 



where 



$(n(z)) = ?io(2)ln(l-?i3(z)) + 

J2niq{z)n2q{z) Y[n2g{z) 



1 - n3(z) 



(l-n3(z))2 



(8) 



describes the excess free energy density due to the 
hard-core interaction. The weighted densities na{z) := 
J2i ^a,i ® gi{z) involve the convolution (0) of the densi- 
ties gi with the weight functions 



i^o,i{z) 


= a{z,Szi), 




Wla;,i(z) 


= SxiCl{z, Szi), 




t^lyAz) 


= Syia{z, Szi), 




i^lzA^) 


^ b{z,Szz), 




'^2x,i{'^) 


— Syib{z, Szi), 




W2j,,i(z) 


— Sxib{z,Szi), 




i^2zAz) 


— SxiSyia{z , Sz 


), 




— SxiSyib{z , Szi 





(9) 

where the abbreviations a(z, S*) :— i((5(-| + z) + (5(-| — z)) 
and b{z,S) := e(| - |z|) are used and Sqi G {0,D} 
denotes the extension of i particles along the q axis. The 
rotational relaxation time 



(10) 
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and the translational diffusion constants 



to /i^ = —1.1 is described in detail. Moreover, the initial 
state discussed here is chosen as 



r 

i^.y 24t' 



367 



(11) 



have been chosen [47|. In the above equations f3 is the 
inverse temperature, A denotes the thermal de Broglie 
wavelength, and rj is the viscosity of the solvent. 

The described DDFT neglects hydrodynamic interac- 
tions between the platelike particles. This is considered 
as a reasonable approximation because of the small par- 
ticle densities in this study close to the bulk isotropic- 
nematic two-phase coexistence region [i^, . 

An equilibrium state fulfills the Euler- Lagrange 
equation 



t^i{z, [g'''^]) = 



(12) 



with the chemical potential /i, i.e., equilibrium density 
profiles render the local chemical potential as a function 
of position (z) and orientation (i) into a constant. There- 
fore the described DDFT is consistent with equilibrium 
DFT because any equilibrium state g°'^ is stationary un- 
der the dynamics represented by Eqs. ©-([S]). Moreover, 
from Eqs. (H]) and one concludes a vanishing particle 
current for equilibrium states: jlg'^'^] = 0. 

For later convenience, the reduced (local) chemical po- 
tential fi* := Pfi — 31og(^) is defined. 



C. Bulk phase behavior and numerical method 

The described model of fluids of platelike colloidal par- 
ticles exhibits a first-order isotropic-nematic bulk phase 
transition at a reduced chemical potential /xj^ = —1.3727. 
At this binodal (b) the isotropic bulk phase of density 
^iso£)3 _ 1.1440 coexists with the nematic bulk phase 
of density g'^°"^D^ = 1.5789 and scalar order parameter 
S^J^"" = (|cos(i9)2 - i) = 0.82696, where (•) denotes the 
thermal average and is the angle between the particle 
orientation and the director 44]. 

Solutions of the DDFT equation specified in Sec. Ill Bl 
are calculated numerically by means of the Euler-forward 
method with integration time steps of 



III. RELAXATION TOWARDS THE STEADY 
STATE 

In this section the relaxation process of an initial state 
g{t = 0) towards the stationary state corresponding to 
the boundary conditions exerted by the two reservoirs 
coupled to the system at z = and z — H (see Sec. |n| is 
studied. Here the arbitrary case of the reservoir at z = 
sustaining the isotropic bulk structure g^ corresponding 
to the reduced chemical potential fj,Q = —1.9 and the 
reservoir at z = H sustaining the nematic bulk structure 
with the director parallel to the z axis g^ corresponding 



e(z,t = 0)= ^« 



,z> 
,z < 



H_ 

2 
H 



(13) 



which approximates the equilibrium structure in the 
channel in the presence of an impermeable membrane 
located at z = -y. Here a detailed discussion of the case 
of the nematic state g^ with the director perpendicular 
to the z axis is not necessary, because the temporal evo- 
lution turns out to be qualitatively the same as for the 
case of parallel alignment. 

The solutions of the DDFT equations Eq. ((2|l- (fTT1) 
have also been calculated for different initial states than 
Eq. p^ . such as the linear density profile g{0 < z < 
H,t = 0) — ^ jj)Qq + jjQjj and the almost empty 
channel p(0 < z < i/, t = 0) « 0, in order to verify that 
the final steady state is independent of the initial state. 

In Fig. [T] the temporal evolution of the initial state 
Eq. (fT3)) for the channel length H = BOD is displayed. 

The total number density g{z,t) is shown in Fig [T^. 
Up to time t ^ lO'^r the initial steplike density pro- 
file smoothens by diffusion until the whole channel at 
z e [0,H] is affected. g{z,t) for times t > lO^r does 
not undergo visible changes, i.e., from then on the steady 
state is (practically) attained. The equilibrium bulk den- 
sity values of the isotropic binodal g'^° and the nematic 
binodal are shown in Fig. [1^, too. The steady state 
is approximately linear below the isotropic and above the 
nematic binodal densities. In the density range between 
the binodals, corresponding to the bulk two-phase coex- 
istence region, a steep portion of the steady state density 
profile is related to the isotropic-nematic bulk phase tran- 
sition. It has been verified by solving the DDFT equa- 
tions that in the case of two reservoirs of both isotropic 
structure a similar steep portion is absent. 

Figure [TJd displays the temporal evolution of the ne- 
matic order parameter tensor component parallel to the z 
axis Qzz{z, t). Similar to the total number density profile 
g{z,t) in Fig.[T^, a diffusive smoothening of the initially 
steplike profile takes place at times t ^ lO'^r whereas 
the steady state has been attained at times t > lO^r. 
Between the order parameter values at the isotropic bin- 
odal and the nematic binodal S^°™ the order parameter 
profile increases rapidly as a function of position. 

For times larger than the rotational relaxation time, 
t ^ r, the local chemical potentials fi* for i parti- 
cles are practically independent of i. Hence it is use- 
ful to consider the reduced chemical potential profile 
fj,*{z,t) := ^Yl,ifJ-*iiz,t), which is depicted in Fig. [1};. 
Again a smoothening of the initial steplike profile in the 
time range t ^ 10** r due to diffusion is followed by a re- 
structuring into the final steady state which is has been 
reached at times t > IO^t. The slope of the steady state 
profile as a function of position z slightly decreases upon 
increasing z. 

In Fig. [TJi the total particle current j(z, t) := 
'^i ji{z,t) is displayed. At early times the current is 
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FIG. 1: (a) Total number density Q{z,t), (b) nematic order paramter tensor component parallel to the z axis Qzz{z,t), (c) 
reduced local chemical potential n*{z,t), and (d) total particle current j{z,t) of platelike colloidal particles of size _D in a 
channel of length H = 50D connecting two particle reservoirs which sustain bulk equilibrium structures corresponding to the 
reduced chemical potentials /xq = —1.9 at z = and (1% — —1.1 at z = H. The structure in the reservoir at z = is isotropic 
whereas the structure in the reservoir at z = H is nematic with the director parallel to the z axis. The isotropic and the 
nematic binodals of the colloidal fluid in equilibrium are displayed for comparison. The time is given in terms of the rotational 
relaxation time r. No further changes are visible for times t > lO'^r. 



localized near the position z = of the discontinuity of 
the initial state Eq. Because /ig < /i^ the current 

is negative. With time the spatial current distribution 
j{z, t) broadens and ultimately becomes homogeneous in 
the steady state. Under the present conditions (reservoir 
chemical potentials, nematic director ahgnment, chan- 
nel length) the steady state current j{z,t = oo)D^t — 
-7.219 • 10"'* is attained. 

In order to quantify the "distance" of a given state g 
from the steady state, which is characterized by a ho- 
mogeneous particle current, consider the nonstationarity 
parameter 



e[Q\ 



max j{z, [g]) 
ze[Q,H] — 



niin 7-(z, [(?]), 

ze[a,H] - 



(14) 



which is a functional of the state g. Obviously e[g] > 0. 
The nonstationarity parameter vanishes if and only if the 
current is homogeneous, e.g., if g is a steady state. 

Figure [2] depicts the nonstationarity parameter e[e(i)] 
evaluated for the state g{t) evolved from the initial state 
Eq. (fT3|) as a function of time t for channel lengths 
H = lOOD (sohd Hue), H = SOD (dashed line), and 
H = 25D (doted line). In Fig. [2^ one identifies a power 
law behavior e ~ i~ 2 within a time range t < t t^ (H), 
where t^ (H) denotes an (approximate) crossover time. 



which depends on the channel length H. The expo- 
nent is a consequence of the diffusive smoothening (see 
Fig. [T|) of the initial steplike chemical potential profile: 
If d{t) denotes the "width" of the interval in Fig. [TJ; 
where ^*{z,t) deviates considerably from the reservoir 
values one expects a diffusive power law d{t) ~ is as 
long as d{t) < H. Since the particle current j{z,t) is de- 
termined by the gradient of the local chemical potential 
(see Eq. (|4])), which in the smoothened interval is approx- 



imately given by 



d{t) 



-, one finds £[e(i)] ^ d{ty 



t 2 . The crossover time t^ (H), which marks the end of 
the diffusive smoothening process, could be defined by 
d{t^{H)) = H, i.e., one expects, t^{H) ~ H"^. At times 
t ^ [H] an ultimate exponential decay of e[£'(i)] is 
found (Fig. [2)3). The decay times t*{H) of this exponen- 
tial decay are t*(100L») = 63260t, t*{50D) = 15793r, 
and T*{25D) = 3959r for the reservoir configurations 
discussed in the present section. These values are pro- 
portional to and they are of the same order of magni- 

tude as — , the one-dimensional translational diffu- 

^t x.y,z 

sion time which is expected to be related to the ultimate 
relaxation time scale. 
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FIG. 2: Nonstationarity parameter £[£i(t)], defined in the main 
text, of state Q{t) at time t for channel lengths H = lOOD 
(solid line), H~= 50D (dashed line), and H = 25D (dotted 
line). The time is given in terms of the rotational relaxation 
time r. In the time range r < t < (H) with t^{H) ~ 
diff'usive smoothening of the initial steplike state leads 
to a power law behavior e ~ (see (a)). At times t ~ 
(H) a crossover to an ultimately exponential decay of e (see 
(b)) takes place. The time scale of this exponential decay is 
proportional to H^. 



IV. PROPERTIES OF THE STEADY STATE 

In the previous section one specific set of reservoir con- 
figurations has been fixed arbitrarily in order to study 
general features of the temporal relaxation towards the 
corresponding stationary state. Here only properties of 
stationary states and their dependence on the reservoir 
configurations is studied. However, the investigation 
is restricted to cases with one reservoir sustaining an 
isotropic structure at z = and one reservoir sustain- 
ing a nematic structure g^^ aX z = H . In contrast to the 
previous section both possible alignments of the nematic 
director with respect to the z axis — parallel and per- 
pendicular — are considered explicitly. We focus on the 
isotropic-nematic interface and on the magnitude of the 
particle current as functions of the boundary conditions 
set by the reservoirs. 

As a general qualitative feature, revealed in the previ- 
ous section, the stationary state varies rapidly with po- 
sition if it is "close" to the (equilibrium) bulk two- phase 
coexistence region (see profiles t — lO^r in Figs. [T^ and 




FIG. 3: Position of the nonequilibrium isotropic-nematic in- 
terface 2'"* in a channel of length H = 50-D with the ne- 
matic director aligned parallel and perpendicular to the z 
axis. In (a) the isotropic reservoir is fixed to {qo — gb°)D^ — 
—0.04 whereas in (b) the nematic reservoir is fixed to [qh — 
qT'")D^ = 0.4. 

[Ud). This feature is reminiscent of a free interface be- 
tween an isotropic and a nematic phase at equilibrium 
bulk coexistence. Here, however, the steplike profiles are 
nonequilibrium structures; we call them nonequilibrium 
interfaces. There are several possibilities to describe the 
"position" z'"* of the nonequilibrium interface. Here the 
following three are considered: 

e(zf,i = oo) i(p'- + g--), 
{inaxQu){ziS\t^oo) 

i ^ 

/.*(z;"',i = oo) /i^ (15) 

The different definitions of the interface position z'"* 
are displayed in Fig. [3] for a channel of length H — 50D 
and for different alignments of the nematic director of the 
reservoir state g^. They are ordered as z™* < z™* < Zg* 
with zi"* - zl" ~« O.ID. Moreover, z^* - z|"' « 0.3D for 
parallel and Zg* — z^"* « 0.61? for perpendicular align- 
ment of the nematic director. These differences compared 
to the channel length are small, however, hence one con- 
cludes that the three definitions for the interface position 
Eq. are equally reasonable. 

Figure [3^ depicts the interface position z'"* as a func- 
tion of the density in the nematic reservoir gn with 
the density in the isotropic reservoir fixed to goD^ = 
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FIG. 4: Stationary state current j(t = oo) in a channel of 
length H = SOD with the nematic director aligned parallel 
and perpendicular to the z axis. In (a) the isotropic reservoir 
density is fixed to {qq — gf°)D^ = —0.04 whereas in (b) the 
nematic reservoir density is fixed to {qh — Qb^^)D'^ — 0.4. The 
main plots exhibit j{t = oo) as a function of (a) qh and (b) 
qq. An alternative representation of j{t = oo) as a function 
of the chemical potential difference of the reservoirs — fig 
is displayed in the insets. 



giso£)3 0.04. Upon increasing the nematic reservoir 
density qh the nonequihbrium interface is shifted to- 
wards the isotropic reservoir {z ~ 0). Conversely, Fig.[3|3 
displays the interface position z'"* as a function of the 
isotropic reservoir density go where the nematic reser- 
voir density is fixed to QhD^ = g^'^'^D^ + 0.4. Under 
these conditions, the nonequihbrium interface shifts to- 
wards the nematic reservoir (z — H) upon decreasing the 
density of the isotropic reservoir qq. Moreover, for given 
reservoir densities, the nonequihbrium interface for par- 
allel alignment of the nematic director is located closer 
to the nematic reservoir (z = H) than for perpendicular 
alignment. 

The dependence of the total current of the stationary 
state j(t = oo) on the configurations of the reservoirs 
is displayed in Fig. [D In Fig. 0^ the isotropic reser- 
voir density is fixed to goD^ = g'^°D^ — 0.04 whereas in 
Fig. the nematic reservoir density is fixed to qhD^ = 
^ncm£)3 _|_ 4_ xhe main plots exhibit j{t = oo) as a 
function of the reservoir densities gu (Fig. [3^) and go 
(Fig. [4)d); the insets alternatively display j{t = oo) as a 
function of the chemical potential difference of the reser- 
voirs iJL*fj — /Xq. Since qh > go, or equivalently /i^ > /ip, 
the total current is negative. The magnitude \j{t = oo)\ 



increases with qh — go and it is larger for perpendicular 
than for parallel alignment of the nematic director. The 
latter observation is not surprising if one recalls T^^y > F^ 
(Eq. pT|) ). However, the two cases shown in FigsS^ and 
|4}d are remarkably different in the sense, that the cur- 
rent difference between perpendicular and parallel align- 
ment \j±(t = oo) \ — \ j\\it = oo)\ is almost constant upon 
fixing the nematic reservoir density gn and varying the 
isotropic reservoir density go (Fig. \^) whereas this dif- 
ference increases upon increasing the nematic reservoir 
density gn and fixing the isotropic reservoir density go 
(Fig. 2^) . Finally, for the same reservoir density differ- 
ence gn — Qo, the magnitude of the stationary state cur- 
rent \j{t = oo)\ in Fig. |3)d where most of the channel is 
filled with isotropic fluid is larger than in Fig. 3^ where 
the channel contains predominantly nematic fluid. 



V. DISCUSSION AND SUMMARY 

In the present work nonequihbrium steady states of a 
fluid of platelike colloidal particles in a channel which 
connects two reservoirs sustaining bulk structures of dif- 
ferent chemical potentials have been investigated. A typ- 
ical platelet fluid is sample AlOP of Ref. which con- 
sists of an aqueous dispersion (solvent viscosity 77 = 8.9 • 
10~^Pa • s) of sterically stabilized gibbsite platelets of di- 
ameter D = 165nm. According to Eq. ^TU\i the rotational 
relaxation time is given by t ~ 214/is. The calculated 
translational diffusion coefficients (seeEq. ([rTjl ) compare 
well with the measured ones of Ref. [50|. The numeri- 
cal results of Sec. IIIII suggest that the fluid in a chan- 
nel of length H — 50D « 8.3/im has definitely reached 
the steady state within a time t = lO^r « 21s. The 
corresponding steady state current is \j{t = oo)\D'^ — 
7.219 • 10"^/r « 3.4s-^ In this work the channel width 
is assumed to be much larger than the particle size such 
that effects of the channel walls can be neglected. For a 
channel width of 10/zm « 61Z?, say, a steady state chan- 
nel current of approximately 12400 platelets per second is 
found. As a reservoir of volume 1£ = O.OOlm'^ contains of 
the order of 10^^ platelike particles, the reservoir density 
in an experimental realization of the setting discussed 
here stays constant over any conceivable experimental 
time scale. 

It has been shown in Sec. IIIII that the relaxation pro- 
cess of the platelet fluid towards the steady state is purely 
diffusive. It comprises a diffusive smoothening of the ini- 
tial steplike fluid structure until the complete channel is 
affected, which takes a time proportional to H^, followed 
by a structural relaxation corresponding to the slowest 
diffusion mode with a relaxation time proportional to 
H^. The two different regimes have been detected by 
a power law and an exponential decay, respectively, of 
the nonstationarity parameter e (Eq. (|14[) ). which mea- 
sures the inhomogeneity of the particle current within the 
channel. In the light of the rather complicated nonlocal 
fluid model (Sec.|Tl| involving translational as well as ori- 
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entational degrees of freedom, the clear identification of 
the power law and exponential decay regimes of e as a 
function of time suggests that these regimes might also 
be easily found in real platelet fluids. 

In summary the present work studies the formation 
and the structure of nonequilibrium steady states in flu- 
ids of platelike colloidal particles in a channel by means 
of dynamic density functional theory. Localized rapid 
changes of density (Fig. [1^) and order parameter tensor 
(Fig. [TJ)) profiles of nonequilibrium steady states are sim- 
ilar to free interfaces. The local chemical potential pro- 
file interpolates smoothly between the values sustained 
by reservoirs at the ends of the channel (Fig. [TJ;). The 
broadening of the particle current distribution (Fig. [TJl) 
suggests the introduction of a parameter measuring the 
distance from the steady state. The purely diffusive re- 
laxation process towards the steady state comprises two 
regimes: a smoothening of the initial steplike structure 



followed by an ultimate relaxation of the slowest diffusive 
mode (Fig. [2]). The position of a nonequilibrium interface 
(Fig. [3]) and the particle current (Fig. [4|) of steady states 
depend nontrivially on the structure of the reservoirs due 
to the coupling between translational and orientational 
degrees of freedom of the fluid. 
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